Finite temperature inelastic mean free path and quasiparticle lifetime in graphene 
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We adopt the GW approximation and random phase approximation to study finite temperature 
effects on the inelastic mean free path and quasiparticle lifetime by directly calculating the imagi- 
nary part of the finite temperature self-energy induced by electron-electron interaction in extrinsic 
and intrinsic graphene. In particular, we provide the density-dependent leading order temperature 
correction to the inelastic scattering rate for both single-layer and double-layer graphene systems. 
We find that the inelastic mean free path is strongly influenced by finite-temperature effects. We 
present the similarity and the difference between graphene with linear chiral band dispersion and 
conventional two dimensional electron systems with parabolic band dispersion. We also compare 
the calculated finite temperature inelastic scattering length with the elastic scattering length due 
to Coulomb disorder, and comment on the prospects for quantum interference effects showing up 
in low density graphene transport. We also carry out inelastic scattering calculation for electron- 
phonon interaction, which by itself gives rather long carrier mean free paths and lifetimes since 
the deformation potential coupling is weak in graphene, and therefore electron-phonon interaction 
contributes significantly to the inelastic scattering only at relatively high temperatures. 



I. INTRODUCTION 

Ever since the pioneering work by Novoselov and 
Geimi, graphene has attracted attention because of 
its potential application in future nanotechnology. Its 
chiral 2D linear dispersion leads to novel proper- 
ties, distinguishing it from conventional parabolic band 
semiconductors. Electron-electron (e-e) interaction in 
graphene is of great interest, both experimentally 2 ^— and 
theoretically^—, because it plays an important role in 
determining electronic properties of graphene. The re- 
duced dimensionality dramatically enhances the electron- 
electron interaction effects^ while chirality and linear 
dispersion have subtle qualitative and quantitative ef- 
fects distinguishing graphene from the usual non-chiral 
parabolic semiconductor-based 2D systems. 

Experimentally, the inelastic mean free path of 
graphene due to e-e interaction, denoted as I throughout 
this paper, is an essential parameter relevant for possible 
ultra-fast device applications because this length defines 
the distance over which the carriers move without any 
energy loss^. In addition, several experimental groups 
have investigated the temperature and gate-voltage de- 
pendence of phase coherence length in graphene^—, 
which is directly determined by I, although may not al- 
ways be identical to it. Thus, I plays an important role 
in determining quantum interference induced localization 
phenomenon. The electro-electron interaction induced 
inelastic scattering is a dominant factor in determin- 
ing the weak localization effects at low temperatures^ 
since the phase coherence length mainly originates from 
Coulomb interaction among electrons. In particular, the 
inelastic mean free path or more specifically whether it 
is longer or shorter than the elastic transport mean free 
path is an important ingredient in understanding the ori- 
gin of the recently observed metal-insulator transition in 
high-purity graphene devices^—. 



Theoretically, the inelastic quasiparticle lifetime of ex- 
trinsic or doped graphene at zero temperature has been 
considered in the literatur e) 20 ' 21 . (In intrinsic undoped 
graphene, the inelastic quasiparticle lifetime due to e- 
e interaction vanishes at the Dirac point, indicating the 
Dirac point to be a non- Fermi-liquid unstable points, but 
this is not a problem at any finite carrier density, which is 
the case we primarily study in this paper.) In doped ex- 
trinsic graphene, it has been found that the chiral linear 
dispersion of graphene results in qualitative differences 
in the energy dependence of inelastic quasiparticle scat- 
tering rates compared with conventional parabolic two 
dimensional electron systems (2DES). Recently, Schiitt 
et at. have analyzed the inelastic scattering rate in- 
duced by Coulomb interaction for intrinsic, i.e. undoped 
graphenoS. They have shown that finite temperature 
strongly affects the inelastic quasiparticle lifetime at the 
Dirac point. The inelastic scattering rate corresponds 
to the so-called on-shell imaginary part of self energy 
ImE fl (k, £k), k being the momentum and £k the quasi- 
particle energy measured from the chemical potential /a. 
So far, neither the inelastic mean free path nor the imagi- 
nary part of the self-energy of extrinsic graphene at finite 
temperature has been explicitly calculated. In addition, 
the effect of finite temperature on the off-shell imaginary 
part of the self-energy of intrinsic graphene has not yet 
been studied. 

In this paper, we theoretically study the inelastic mean 
free path I and the imaginary part of the retarded self- 
energy ImS fl (k, w) in graphene within the leading or- 
der many-body perturbative GW approximation. The 
GW approximation involves the leading-order theory in 
the dynamically screened Coulomb interaction which in- 
cludes all the ring diagrams in the electron self-energy. 
We generalize the previous GW work2i in the literature 
(carried out at T = 0) to finite temperatures, for both 
extrinsic and intrinsic graphene. We find that the finite 
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temperature results for I and ImS fl (k, lo) are very dif- 
ferent from the zero temperature ones due to a dramatic 
change in the dynamic dielectric function with increasing 
temperature. We also provide an analytical expression 
for the on-shell imaginary part of the self-energy for a 
double-layer graphene system, which is relevant for a re- 
cent experiment^. We also calculate the imaginary part 
of the electron self-energy due to electron-phonon (e-p) 
interaction in graphene through the deformation poten- 
tial coupling, and provide results for the corresponding 
inelastic mean free path as a function of temperature. 
In general, phonon effects are rather small in graphene 
because the deformation potential coupling is very weak 
in graphene. Thus, except at very high temperatures, 
phonon effects are negligible in the determination of the 
inelastic mean free path and quasiparticle lifetime. 

The rest of this paper is organized as follows. In Scc.lnl 
we introduce the theoretical formalism for treating inelas- 
tic e-e scattering at finite temperatures in both mono- 
layer and double-layer graphene systems. We also obtain 
analytical expressions for the imaginary part of the re- 
tarded self-energy in the low energy and low temperature 
limit. In Sec. lIIII we provide numerical calculations of the 
imaginary part of the self-energy (which is the inverse 
of the quasiparticle lifetime) and the associated inelas- 
tic mean free path as a function of carrier energy and 
temperature. We also compare the inelastic mean free 
path with the disorder induced elastic mean free path. 
In Sec. IIV1 we provide the electron-phonon interaction 
results separately, giving both the basic theory briefly 
and the numerical results. In Sec. [Vj we discuss and 
summarize the main results of this paper. 



the Coulomb interaction in momentum space with aver- 
age background dielectric constant n. s, s' = ±1 are band 
indices. £k+q,s = s£k+q — (£k = ft«F|k| with graphene 
Fermi velocity Vp = 10 6 m/s) is the carrier energy rel- 
ative to the finite temperature, non-interacting chemical 
potential /i, determined by the conservation of the to- 
tal electron density. e(q,u>) = 1 + V q IL(q, lo) is the finite 
temperature dynamic dielectric function within random- 
phase approximation (RPA) 22 i 23 , and II(g, lo) is the irre- 
ducible polarizability. Note that both single particle ex- 
citations (Im[e] 7^ 0) and plasmon excitations (Re[e] = 
and |Im[e]| = + ) contribute to ImSf (k,uj). We use 
H = 1 throughout. 

The inverse quasiparticle lifetime or the quasiparticle 
scattering rate 1/r = T s (k) is directly related to the 
imaginary part of the on-shell self-energy, i.e., 

i = r s (k) = 2imsf(k,a, s ) (2) 

T 

The inelastic mean free path I due to e-e interaction is 
then given by 

I = VpT (3) 

Thus the evaluation of Im£f (k, £k,s) immediately leads 
to the calculation of inelastic mean free path and quasi- 
particle lifetime. We calculate the imaginary part of the 
self-energy in details in this paper. 



II. THEORETICAL FORMALISM 
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In this section, we provide the theoretical formalism for 
evaluating the finite temperature imaginary part of the 
retarded self-energy due to e-e interaction in graphene. 
We also obtain the asymptotic behavior of the imaginary 
part of the self-energy in the low energy and low temper- 
ature limit. 



A. Imaginary part of self-energy for extrinsic 
doped graphene 

The imaginary part of the retarded self-energy of 
monolayer graphene within GW approximation can be 
expressed as^^: 

ImEf (k,Lo) = -- E v qi n B(£,k+q,s' - w) + rii?(£ k+q , s <)] 

q,s'=± 

(1 + ss'cos#)Im 



1 



efeCk+q,. 



(1) 



where functions np and jig denote the Fermi and Bose 
distribution functions, respectively V q = 2ire 1 / (nq) is 



We first give the analytical formula for ImS^(fc, at 
zero temperature and for ep ^ This has already been 
given in Ref. [f| and we just provide the results below for 
completeness and comparison. 

The quasiparticle scattering rate can be obtained 
through the imaginary part of the self-energy within the 
on-shell approximation lo = £k,s: 

Im£«(fc,60= \ E E - sWq) - *(-&+,,.')] 

s'=± q 

x(l + s'cos0)lm — — (4) 

V ' U((7,s'£k + q-£k+*0+)J V ' 

where i?(x) is the Heaviside unit step function. 

From Eq. <(4j) , it is straightforward to see that the inter- 
band scattering from the valence band s' = — 1 vanishes 
at zero temperature. In addition, the plasmon emission 
process also vanishes in calculating the on-shell imag- 
inary part of the self-energy because the plasmon fre- 
quency requires Lo p i{q) > hvpq for extrinsic graphene at 
zero temperature 2 ^ but |£k+q — £k| < fivpq in the dy- 
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namic dielectric function. Thus, Eq. (j4]) reduces to 

2 p2n rk 

Im£?(/c,£ fc ) = / d9(l + cos9) dk'k' 

47rK Jo Jk F 



;; Im[l/6(|k-k'|,e fc . ~e k +iO+)} 
\/k 2 + fc' 2 - 2kk' cos 9 



In the long- wavelength x — q/2kp <C 1 limit, we can 
obtain the analytical expression for Eq. In this limit, 
the dominant contribution to Im[l/e] comes from low en- 
ergies u — u/vpq <C 1, where the irreducible polarizabil- 
ity in the leading order is given by^ 2 - 



n(g,t 



Dofl 



(6) 



where _D 



o 



2fc F 



, the density of states at Fermi energy. 



Then, the asymptotic behavior of the imaginary part 
of the dynamic dielectric function can be written as: 



Im- 
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e(q, e k +q - £k + iCH 



fc' - fc 
qTF 



(7) 



where fc' = |k+q| and qxF — 4 frf^ K denotes the Thomas- 
Fermi screening wave- vector in graphene. To obtain 
Eq. (|7J|, we have used the fact that the integrand in 
Eq. <j5j) has sharp peaks near forward scattering momen- 
tum transfer, i.e. q = 0. The asymptotic expansion of 
the imaginary part of the self-energy in Eq. (|S|) can then 
be obtained by integrating out the angular part and ex- 
panding in the small momentum transfer fc — fc' fc^: 
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(8) 



The leading order energy correction to the on-shell imag- 
inary part of the self-energy is £f ln£fc, similar to that 
obtained for conventional two-dimensional electron sys- 
tems (2DES)^. The subleading term in Eq. © has op- 
posite sign to that in 2DES because of the chiral nature 
of graphene and the absence of interband contribution 
from valence band (s' = — 1) at zero temperature. 



2. Finite temperature analytical results T > 



To maintain analytic tractability for Im£^ + at finite 
temperature, we consider only low temperatures and 
small i.e., Ef S> k B T 3> £fc- In this limit, both the 
interband scattering and the temperature dependence of 
Im[l/e], being exponentially suppressed at low tempera- 
tures, can be ignored. Applying the same technique as in 
the T = case, the finite temperature on-shell imaginary 



part of the self-energy can be approximated as: 

,2 c2tt 



e' 

47TK 

n F (sk)] 



d6(l + cosO) / dk'k'[n B {e k , - e k ) 
Jo 

Im[l/e(|k-k'|,e fc / -e k + iQ+)} 



Vfc 2 + fc' 2 - 2fcfc' cos 9 



(9) 



The imaginary part of the inverse dielectric function can 
still be approximated as in Eq. (|7|). After integrating out 
the angular part and expanding in the limit £ k <C £f, 
the momentum integration is very straightforward and 
we obtain an asymptotic formula for Eq. (|9]) in the limit 
e F > k B T^£ k : 



ImS 



k+ 



AksTf 
8s F 



W( kBT ) 
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(10) 



Thus, the leading order temperature correction to the 
on-shell ImS^ + goes as t 2 \nt with t = T/T F (T F being 
the Fermi temperature). 



B. Imaginary part of self-energy for double- layer 
graphene 

In this subsection, we provide analytical results for 
ImSjJ, of a double-layer graphene system consisting of a 
"studied" and a "control" layer of graphene separated by 
an insulating layer—, which has recently been studied ex- 
perimentally. These two layers of graphene are electroni- 
cally isolated but the screening effect of the control layer 
must be considered. We will only focus on the doped 
double-layer case, where Fermi levels of both layers are 
above the Dirac point. 

The imaginary part of the self-energy of the "stud- 
ied" layer for a double-layer graphene, in the presence of 
screening by carriers both in the control layer and the 
studied layer itself, can be written as 2 ^: 



ImEf(fc, £ 



'=±.q 



k+q,s' 



k+q,s' 



x (1 + ss' cos 0)lm 



k+q,s' 



)}V(q) 



(11) 



e scr (q, £k+q.s' — w) is the RPA dynamic dielectric function 
for the double-layer system, which incorporates screening 
effects from a nearby "control" graphene layer in addi- 
tion to the usual screening by the studied layer electrons 
themselves 2 ^. 



1 



l + V(q)IL 22 {l-e- 2 i d ) 



e s -(g, u) 1 + V(q)[U n + n 22 ] + n u Tl 22 V(q) 2 [l - e' 2 ' 

(12) 

where IIn and 1I 2 2 are polarization operators of the stud- 
ied and the control layer, respectively, d is the separation 
between the two layers. Note that Eq. (fT2)) reduces to the 
dynamic dielectric function of monolayer graphene in the 
large d limit as it should because the two layers become 
decoupled in the large d limit. 
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1. Zero temperature analytical results T = 

In this section, we derive the analytical formula for 
the imaginary part of the double-layer self-energy at zero 
temperature in the limit sp ^> £fc. This is very close to 
what we derived in the previous subsection for monolayer 
graphene but with different prefactors coming from the 
additional screening effects of the nearby second layer. 

At zero temperature, the interband scattering from the 
valence band s' = — 1 vanishes in Eq. (fTTj). and the on- 
shell expression reduces to : 



Im£*_ 



2jt 



/ d9(l 

4ttk Jo 



cos ( 



dk'k' 



Im[l/6 flcr (|k-k'|,e fc , -e fc +iO+)] 
Vk 2 + k' 2 - 2kk' cos 9 



(13) 



In the long- wavelength x = q/2kpn i 2) <C 1 and low 
energy limit u = u/vpq <C 1, it is possible to perform an 
analytical evaluation of Eq. (fT3")) . In this case, the irre- 
ducible polarizability of graphene is given in the leading 
order by^ 



IIu^, w) — Am[l + i ] 

v F q 

UJ 

n 22 ((?,w) ~ L»022 1 + i 

VFq 



(14) 
(15) 



where -Don = and D 022 = correspond to the 

density of states at Fermi energy of the studied and the 
control layer, respectively. 

Then, the imaginary part of the screened dielectric 
function is approximated as: 



Im 



1 



-(jfe'-fc) 



e scr (q,e k+(l -e k +iO+) 

qTFl + qTF2 + ^qTFiqTF2<i + 4q T Fl{qTF2d) 2 



where qrFi 



(qTFi 

_ A k F1 e 



9TF2 + 2qTFiqTF2d)' 2 



(16) 



and qrF2 



4A- 



E2S - are the 



Thomas Fermi wave-vectors of the studied and the con- 



trol graphene layer. In deriving Eq. (fTr?)) . we have used 
the long- wave length limit |k — k'|/2fc_p <C 1, which also 
indicates that |k — k'|d <C 1 for d ~ 0(1 nm) and 

<C 1. From Eq. (fl"6|) . we see that the difference 
in the analytical expressions for ImS^ between mono- 
layer and double-layer graphene comes from the second 
line of Eq. ([T6|l . 

Using the approximation of Eq. (|16[) in Eq. fp~3|) . we 
obtain the final analytical formula: 




(1 + y/ri^/ ^fn{ + 8r s d^ 



(17) 



where n% , n 2 are the carrier densities of the studied layer 
and the control layer respectively. r s = e /(Jvufk) is the 



graphene fine structure constant defining the electron- 
electron interaction strength and Efi is the Fermi en- 
ergy of the studied layer. It is clear that Eq. (fTT)) , as 
expected, reduces to the monolayer graphene result as 
given in Eq. (|SJ) for two extremes: n 2 = and d — > 00. 



2. Finite temperature analytical results T > 

Next, we provide the analytical expression for the 
imaginary part of the self-energy of a double-layer 
graphene in the low T and small £ limit, i.e., £f(i,2) ^ 
ksT >• £fc as before. 

In the low temperature limit, we neglect both the in- 
terband scattering and the temperature-dependent imag- 
inary part of the dielectric function, which are exponen- 
tially suppressed. Then, the asymptotic expansion of 
Eq. (fTTj) is given by: 



ImEf + 



7r(fc B T) 2 



In 



8efi 



k B T 



1.08387 



8e.fi . 
im 2 + 64:r 2 d 2 irn 2 



(1 + y/na/ y/n~[ + 8r s d^7m 2 )' 



(18) 



C. Imaginary part of self-energy for intrinsic 
graphene 

The imaginary part of the self-energy for undoped or 
intrinsic monolayer graphene can be obtained by setting 
fj, = in Eq. (JTJ), in which case the chemical potential is 
always independent of temperature. 

At zero temperature, the analytical form of the irre- 
ducible polarizability for intrinsic graphene^ is: 



n = <f_ I 0(v F q-oj) 
4 \y/(vj,q 2 -u 2 ) 




(19) 



Then, the imaginary part of the self-energy within the on- 
shell approximation s = at zero temperature reduces 
to the product of two # functions^: 

ImS«(fc, e k ) ~ £ q tf(|e k+q - e k | - e q )d(e k - e k+q ) (20) 

The first $ function is due to the imaginary part of the 
dynamic dielectric function while the second one comes 
from the sum of the Fermi and Bose distribution func- 
tions at T = 0. It is easy to see from the above equation 
that ImS^fc, Ek) vanishes because of phase space restric- 
tions imposed by the 1? function. On the other hand, we 
can evaluate the zero temperature imaginary part of the 
self-energy at zero momentum k — and finite energy ui 
by using the polarizability given in Eq. (|19l) . 

Then the imaginary part of the self-energy at k = is 
given by: 



ImE£(fc = 0,w)=w/(r a ) 



(21) 
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with 



/(*) 



2 / f . , , 2-(irx/2)\ nx-^rx 

—5— nil — x) H = m == — 

7v x V V ( 7ra; ) 2 _ 4 ^ + V C 7 ™) 2 _ 4 

n (22) 

which was already obtained in Ref. [g . Note that for the 
general case (k, w), it is almost impossible to get an an- 
alytical formula for the imaginary part of the self-energy 
due to the complicated angular integration. However, we 
know that ImE5(k, u>) = for u> < hvpk at zero temper- 
ature because of the combination of the two •& functions 
mentioned above. 

Next, we mention that the analytical formula for the fi- 
nite temperature imaginary part of the self-energy for the 
intrinsic undoped graphene is quite tricky necessitating 
a very careful analysis of the finite temperature dynamic 
dielectric function. Since undoped graphene is not the 
focus of our current work, we refer the reader to Ref. Q , 
where the on-shell imaginary part of the self-energy in 
intrinsic graphene has been recently studied. Note that 
the finite temperature imaginary part of the on-shell self- 
energy is proportional to the temperature T, as can be 
shown by dimensional counting and can also be seen from 
the numerical results shown in Sec. IIIII In addition, we 
find that the leading order temperature correction to the 
imaginary part of the self-energy for two cases: k = 
and ksT <C |cj — £k|, are exponentially suppressed. Since 
the self-energy for the undoped case has already been 
studied analytically in some details in Ref. [9[ we do not 
provide any further discussion of this issue in this paper. 
We also add that from the experimental perspective, pure 
intrinsic graphene is uninteresting since it is unstable to 
any (for example, disorder-induced) density fluctuations 
which would locally dope the Dirac point, and thus, all 
experimentally studied graphene samples are likely to be 
extrinsic doped graphene. 
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III. NUMERICAL RESULTS 

In this section, we provide our numerical results for 
both inelastic mean free path and the imaginary part 
of the self-energy (which is essentially the inelastic scat- 
tering rate or the inverse of the quasiparticle lifetime) 
in graphene for both extrinsic graphene and intrinsic 
graphene. If not specified, the effective background di- 
electric constant used in our calculation is k = 5, i.e. 
r s = 0.44, corresponding to monolayer graphene sand- 
wiched between two Boron- Nitride layers^. We compare 
our numerical results with the asymptotic analytical re- 
sults derived in Sec. HI1 We also calculate the ratio be- 
tween the inelastic mean free path and the elastic mean 
free path induced by long-range Coulomb disorder as a 
function of temperature in order to assess the importance 
of possible quantum interference effects in highly resistive 
graphene samples at low temperatures. 
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FIG. 1. (Color online), (a). Calculated ratio A, which is 
the numerical result of ImE+(kF, £kp) divided by the analyt- 
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(given in Eq. (Tl0|) ), as function of T/Tf- The inset 
presents A in the low temperature regime, (b) Calculated 
ratio a, the numerical result of ImE+(k,£k) divided by 



hi 



8e F 



+ 



1 



given in Eq. 1181 ) as a function 

of (,k/EF at T = 0. The inset presents a in the low en- 
ergy regime. The solid and dashed lines correspond to k = 5 
(r s = 0.44) and k = 1 (r s = 2.2), respectively. 



A. Inelastic mean free path and imaginary part of 
self-energy for extrinsic graphene 

In Fig. [TJ we compare our calculated analytical results 
of the on-shell imaginary part of the self-energy with the 
numerical results, which show good agreement with each 
other. The ratio A is ImE^kp, ^kj?), calculated numeri- 
cally at finite temperature, divided by its value given by 
Eq. (fTU|) . while the ratio a is Im£+(k, £k), calculated nu- 
merically at zero temperature, divided by its value given 
by Eq. ©. Both ratios A and a are closer to unity for 
smaller background dielectric constant k. This is because 
the approximation used in Eq. (|7|) is more accurate for a 
larger value of qTF- Note that both ratios A and a are in- 
dependent of Fermi-energy Ep for a fixed fine structure 
constant r s within the GW approximation. We want 
to mention that both ratios in graphene have opposite 
trends compared with 2DES in the sense that in the lat- 
ter system A and a are both larger than unity as they 
approach unity asymptotically^ whereas in graphene the 
ratios always lie below unity. We mention that it is quite 
interesting to note that the asymptotic expressions, al- 
though they are derived for energy or temperature being 
much smaller than Ep or Tp, appear to be valid quanti- 
tatively within a factor of two well outside the asymptotic 
regime, i.e. even when T ~ Tp. 

Fig. [2] shows our numerical results for the energy and 
temperature dependence of the inelastic mean free path 
I and the associated on-shell imaginary part of the self- 
energy. In Fig. [U[a), we present the inelastic mean free 
path I as a function of temperature for two different car- 
rier energy values. We can see that / is a monotonically 
decreasing function of temperature. In particular, the 
injected electron can only decay via intraband processes 
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at zero temperature as discussed in Sec. Ill A II In addi- 
tion, the decay via plasmon emission is inhibited at zero 
temperature due to phase-space restrictions^. On the 
other hand, at finite temperatures, the injected electron 
can decay via both interband and intraband excitations 
and the region of single particle excitations also increases 
due to thermal smearing effects. Fig. EJb) shows the en- 
ergy dependence of I for different temperatures and car- 
rier densities. At low temperatures, the mean free path of 
hot-electrons (e^ > Ep) is shorter than that of the quasi- 
particle in the vicinity of the Fermi energy. However, 
very interestingly, the energy dependence of I becomes 
non-monotonic at higher temperatures (see the lowest 
dashed line in Fig. EJb)). This arises from the compe- 
tition between damping through valence band and con- 
duction band, which have opposite trends as a function 
of carrier energy. To be more specific, the scattering rate 
from the valence band contribution (s' = —1) is a mono- 
tonically decreasing function of ^/Ep while the contri- 
bution from the conduction band (s' = 1) is a mono- 
tonically increasing function of ^/Ep. Figs. EJc) and 
(d) provide the on-shell ImS^(k, £k)- We point out that 
ImE£(k, &)/E F for the two carrier densities at zero tem- 
perature in Fig. 0[d) are essentially the same because the 
rescaled ImE^(k/fci? , £k/Ep)/Ep is universal at T = 0. 
Notice that in our theoretical GW formalism, we exclude 
multiparticle excitations, which have a relatively small 
effect on / in graphene^i because of the small values of 
r s implying weak e-e interaction strength. 

For the purpose of comparison, in Fig. [3l we show re- 
sults of / and ImS + for suspended graphene with the 
background dielectric constant k = 1 (r s = 2.2). From 
Figs. [2] and we can see that the results of the mean free 
path I for different values of r s are qualitatively very sim- 
ilar. The inelastic scattering length / is larger for smaller 
r s because of weaker electron-electron interaction. 

To capture the effects of plasmons and interband 
single-particle excitations, we have to consider the off- 
shell imaginary part of the self-energy ImS fi (k, to) (uj ^ 
£k.s), which is also needed to explain the experimen- 
tal ARPES datai^L^i. The numerical results for 
ImE fl (k, u>) are shown in Fig. Unlike the calculation 
for the on-shell imaginary part of the self energy, we can 
see that the off-shell imaginary part of the self-energy 
contains two kinds of contributions from Eq. The 
first comes from the single particle excitation occurring 
at Ime ^ 0. This can be further divided into intraband 
and interband excitations because of the gaplessness of 
graphene. The second contribution comes from plasmon 
excitations, where both Ime and Re e vanish. The plas- 
mon resonance consists of a 5 function at zero tempera- 
ture, of which the width will be broadened as tempera- 
ture increases. Figs. 2] (a) and (b) present the rescaled 
conduction band ImS^ (k, u)/Ep with k = 1.5kp as a 
function of rescaled temperature T/Tp and energy ui/Ep. 
We point out that ImS^(k /kp,ui/Ep)/Ep is universal, 
not depending on the position of the Fermi level but only 
on the graphene fine structure constant r s . The non- 




FIG. 2. (Color online). Calculated inelastic scattering length 
I and the associated on-shell imaginary part of the self-energy 
ImE+ (k, £k)/EF for extrinsic graphene (ef > 0) with dielec- 
tric constant n = 5 (r s = 0.44). (a) I as a function of tem- 
perature T. The solid, dashed and dot-dashed lines are for 
carrier density n = 10 10 , 10 11 and 10 12 cm' 2 , respectively. 
The upper three lines are for = £j, F while the lower three 
lines are for ffc = 0.5ef. (b) I as a function of ^/Ef for differ- 
ent temperatures. The solid and dashed lines are for carrier 
density n = 10 11 , and 10 10 cm -2 , respectively. For the same 
carrier density, the lines from top to bottom correspond to 
T — K, 50 K, and 100 K. (c) and (d) are the associated 
ImE^(fc, ^iJ/Ef corresponding to (a) and (b), respectively. 



monotonic temperature dependence of ImE^: for fixed 
oj/Ep as shown in Fig.[4ja) comes from the competition 
among interband, intraband and plasmon excitations. In 
the zero temperature limit, it has been shown in Ref. pH | 
that the plasmon contribution dominates in the low en- 
ergy regime for k/kp > 1 and it vanishes in the higher 
energy regime. As the temperature goes up, the plas- 
mon contribution is not a 8 function anymore and it gets 
broadened by Landau damping (see Fig. [4] (b)). Also, 
the boundaries of both inter- and intra-band excitations 
are thermally smeared. The above two effects give rise 
to the appearance of a smaller bump in ImS^ in the 
lower energy region at finite temperature as shown in 
Fig. [H (b). Figs. @Jc) and (d) show the imaginary part 
of the self-energy for the valence band ImS^. ImS^ is 
generally an increasing function of temperature. In the 
k = and T = limit, only single particle excitations 
contribute to ImS^(0,w) and there is no plasmon emis- 
sion because it requires q < kp. But the integrand in 
Eq. ([T]) has non-zero values only for q > kp, being re- 
stricted by the sum of the Bose and Fermi distribution 
functions, i.e., (uj + Ep — e q ){Ep — e q ) < 0. In particular, 
only intraband single particle excitations contribute for 
energy ujjEp < 1 and the interband single particle exci- 
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FIG. 3. (Color online). Calculated inelastic scattering length 
I and the associated on-shell imaginary part of the self-energy 
ImE^(fe, £k)/EF for extrinsic graphene (ef > 0) with dielec- 
tric constant n = 1 (r s = 2.2). (a) I as a function of tempera- 
ture T. The solid, dashed and dot-dashed lines are for carrier 
density n = 10 10 , 10 11 and 10 12 cm -2 , respectively. The up- 
per three lines are for £k = £fc F while the lower three lines 
are for = 0.5eF. (b) I as a function of ^/Ef for differ- 
ent temperatures. The solid and dashed lines are for carrier 
density n = 10 11 , and 10 10 cm -2 , respectively. For the same 
carrier density, the lines from top to bottom correspond to 
T — K, 50 K, and 100 K. (c) and (d) are the associated 
ImE+ (fe, corresponding to (a) and (b), respectively. 




T/T p (0/E F 



FIG. 4. (Color online). Calculated ImS R (k,u}) for extrinsic 
graphene (the chemical potential [i > 0) with k = 5. (a) and 

(b) correspond to conduction band while (c) and (d) corre- 
spond to valence band. In Fig. (a), the solid, dashed and 
dot-dashed lines correspond to uj/Ef = 0.5, 1.0 and 2.0, re- 
spectively, (b) The solid, dashed and dot-dashed lines cor- 
respond to T/Tf = 0, 0.2 and 0.4, respectively. In Figs. 

(c) and (d), the solid and dashed lines correspond to k = 
and k = 1.5fc_F. (c) For the same momentum, the lines from 
bottom to top correspond to uj/Ef = 0.5, 1.0 and 2.0, respec- 
tively, (d) For the same momentum, the lines from bottom to 
top correspond to T/Tf = 0, 0.2 and 0.4, respectively. Note 
that ImE^fc = 0,w) = lmE R (k = 0,w). 



tation contribution increases sharply around oj/Ep ~ 1. 
On the other hand, for k = but at nonzero tempera- 
ture, lm£^(0, oj) becomes much smoother compared with 
the zero temperature case. In addition, ImE^OjO) has 
a finite value increasing with temperatures due to the 
thermally smeared boundaries of single particle excita- 
tion continua. Note that ImS 77 = ImE^ for k = as 
seen from Eq. (JXJ) , in which case the equation has no 
band- index dependence. 



B. Comparison between inelastic and elastic mean 
free path in extrinsic graphene 

In this subsection, we compare the inelastic mean free 
path I induced by e-e interaction with the elastic mean 
free path l e due to charged impurity in the environment. 
The ratio £ = l/l e 1 is the necessary condition for 
quantum interference effects to be operational in exper- 
iments since interference necessarily requires the phase 
coherence of energy eigenstates. A detailed discussion 
on this issue has been given in Ref. [18j . where, how- 
ever, the inelastic mean free path was not calculated. In 
Fig. [51 we show our calculated ratio ( as a function of 
temperature for different values of the potential fluctua- 



tion (i.e. puddle) parameter s. This parameter is defined 
as the standard deviation of the probability distribution 
of the disorder potential at a given point in the graphene 
plane. It can be tuned using screening by the second 
nearby graphene layer—. A detailed calculation for the 
elastic mean free path l e has been given in Ref. [28| . Note 
that we do not take into account electron-phonon scat- 
tering in the calculation of the mean free path here since 
the temperature (T < 100K) is relatively low and the 
inelastic mean free path can be even shorter if we in- 
clude electron-phonon scattering mechanism. We defer 
our discussion of electron-phonon interaction to section 
IIVI of the paper. 



C. Inelastic mean free path and imaginary part of 
self-energy for intrinsic graphene 

In this subsection, we show our numerical results for 
the inelastic mean free path I and the imaginary part of 
the self-energy for intrinsic graphene. The behavior of in- 
elastic e-e scattering of intrinsic graphene is quite differ- 
ent from that of extrinsic graphene. Intrinsic graphene is 
a marginal Fermi liquid^. By contrast, extrinsic graphene 
is a well-defined and relatively weak-coupling Fermi liq- 
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10.0 



^ 1.0 



0. 




— s = 1 meV 

- s = 5 meV 

— s = 15 meV 

— s = 25 meV 

50 
T(K) 



FIG. 5. (Color online). Calculated ratio £ = l/l e , the inelastic 
mean free path divided by the elastic mean free path, as a 
function of temperature for carrier density n — 10 9 cm" 2 , 
k = 5, and charged impurity density 9 x 10 10 cm" 2 . The solid, 
dotted, dashed and dot-dashed lines correspond to s = 1.0, 
5.0, 15.0 and 25.0 meV, respectively, indicating increasing 
effects of Coulomb disorder induced inhomogeneous puddles. 
The horizontal line indicates where (above the line) quantum 
interference could play a role. 



uid since typically r s < 1 and is independent of carrier 
density. 

In order to compare with the results provided for ex- 
trinsic graphene, we show the finite- T inelastic scatter- 
ing mean free path I for intrinsic graphene in Fig. |Bfa) , 
(b) and the corresponding imaginary part of the self- 
energy in Fig. [(^c) , (d) using similar parameters as used 
for our extrinsic graphene results. Similar to extrin- 
sic graphene, the inelastic mean free path of intrinsic 
graphene is a monotonically decreasing function of tem- 
perature. However, the leading order temperature depen- 
dence of ImE^(k, £k) for intrinsic graphene is linear in 
T as shown in Fig. [51 At zero temperature, ImE^k, £k) 
vanishes due to phase-space restrictions leading to in- 
finite I. On the other hand, finite temperature effects 
give rise to a nonmonotonic, energy dependent I and 
ImS^(k,^k)- IniS+(k,£k) (or the quasiparticle scatter- 
ing rate) first increases linearly with £k and then de- 
creases with £k- This nonmonotonic energy dependent 
feature of Im£+(k, £k) has also been found in Ref. @, 
and can be explained by the peculiar property of the dy- 
namic dielectric function of intrinsic graphene at finite 
temperatures. 

The numerical results for the off-shell lmJ^ R (k, uj) for 
intrinsic graphene are shown in Fig. [7] Fig. [3a) presents 
ImS^(0, uj) as a function of temperature T . Quite differ- 
ent from extrinsic graphene, where the plasmon peaks are 
broadened by Landau damping due to finite temperature 
effects, intrinsic graphene has no plasmon excitations at 
zero temperature as seen from Eq. (|19[) . Note, however, 
that the plasmon modes are thermally restored at finite 
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FIG. 6. (Color online). Calculated inelastic scattering length 
I and the associated on-shell imaginary part of the self-energy 
ImE+(fc,£fc) for intrinsic graphene (ef = 0) with dielectric 
constant k = 5. (a) / as a function of temperature T for 
different energy e^. The solid, dotted, dashed and dot-dashed 
lines are for energy £k = 1 meV, 5 meV, 10 meV, and 20 meV. 
(b) I as a function of energy u for different temperatures. The 
solid, dashed and dot-dashed lines are for temperature T = 10 
K, 50 K, and 100 K. (c) and (d) are the associated ImE+(fe, tk) 
corresponding to (a) and (b), respectively. 



temperatures in intrinsic graphene. We find that there is 
a peak in lm£^(0, uj) as a function of temperature, which 
arises from the appearance of thermal plasmon modes. In 
particular, the position of these peaks is approximately 
proportional to the energy u because the rescaled func- 
tion ImE R (k, u>) is universal with respect to the rescaling, 
t R = Z R /{k B T), k = hv F k/(k B T) and Q = uj/(k B T). 
Only the line of lm£ fi (0,u;) with lj = 1 meV in Fig.^a) 
shows a peak since the position of peaks for higher energy 
lie outside of the temperature regime we are interested 
in. Fig.[7Jb) shows the energy dependence of ImS^(0, u>) 
for several values of temperature. At zero temperature, 
ImE^O, w) is linearly increasing as a function of ui indi- 
cated in Eq. (|21[) . In the lower energy regime, ImE^(0, uj) 
shows a peak associated with the resonance of plasmon 
excitations. Specifically, both the position and the height 
of the peak are roughly proportional to temperature T as 
mentioned above. The finite temperature ImS^(0, uj) ap- 
proaches the zero temperature ImS?(0, uj) as the energy 
uj further increases. According to Fig.[7]Jc), the tempera- 
ture dependence of off-shell conduction band ImS^(fc, uj) 
is exponentially suppressed while the on-shell ImE^(fc, uj) 
is proportional to T at low temperatures. We plot the 
conduction band lmS R (k,uj) as a function of energy in 
Fig. EJd). One interesting feature is that ImE R (k,uj) 
vanishes for uj < £k and it increases rapidly with uj for 
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T (K) (0 (meV) 



FIG. 7. (Color online). Calculated imaginary part of the 
self-energy lmT 1 R (k y uj) for intrinsic graphene (the chemical 
potential fi = 0) with k = 5. (a) and (b), Conduction band 
as a function of temperature T and energy 
to, respectively, (c) and (d), Conduction band ImE+(fc, oj) 
with k — 10 meV/hvF as a function of temperature T and 
energy u), respectively, (e) and (f) Valence band ImE5(fc,a>) 
with k — 10 meV/hvF as a function of temperature T and 
energy u, respectively. The solid, dotted, dashed and dot- 
dashed lines in Figs, (a), (c) and (e) are for energy u = 1 
meV, 5 meV, 10 meV, and 20 meV. While the solid, dotted, 
dashed and dot-dashed lines in Figs, (b), (d) and (f) are for 
temperature T = K, 10 K, 50 K, and 100 K. Note that 
ImE5(fc = 0,w) = ImE?(fc = 0,u>). 



lj > £k at zero temperature. As the temperature in- 
creases, two characteristic peaks show up in Im£+(fc,w) 
similar to extrinsic graphene as shown in Fig. HJb). This 
arises from an interplay between plasmon enhancement 
and thermally smeared single particle continuum bound- 
aries. We have found that the positions of these peaks 
\<-<Jpi — £k| ar c proportional to the temperature T. We 
provide the valence band Im£^(fc, w) as a function of 
temperature T and energy lj in Figs.[7][e) and (f), respec- 
tively. The low temperature correction to Im£^(fc,cj) is 
exponentially suppressed for uj > due to the constraint 
imposed by the Fermi and Bosc distribution function in 
Eq. ([T]). lmS^(k,oj) increases rapidly at higher temper- 
atures, especially for small energy lj. From Fig. 



we can see that Im£^(fc, w) increases slowly with lj for 
w > £k at low temperatures and Im£^(fc,cj) shows non- 
monotonic dependence on lj at higher temperature. This 
non-monotonicity can be understood through the com- 
petition between a decrease in valence band contribution 
(s' = — 1) and an increase in conduction band contribu- 
tion (s' = 1) with increasing energy u>. 



IV. MEAN FREE PATH DUE TO 
ELECTRON-PHONON INTERACTION 

In this section, we present results for the mean free 
path l ep induced by electron-phonon interaction. We 
analytically calculate its high-(T 3> Tbg) an d low- 
temperature (T <C Tbg) limits, while in the interme- 
diate temperature regime we provide numerical results. 
In addition, we consider the energy and carrier density 
dependence of l ep for different values of temperature. In 
this work, we consider acoustic phonons only, and neglect 
the optical phonons because their energy in graphene is 
too high (~ 2000 K) for them to play any important role 
in the temperature range (T < 500 K) of our interest. 



A. Imaginary part of self-energy due to 
electron-phonon interaction 

The mean free path l ep is directly related to the 
on-shell imaginary part of the self-energy, i.e. l ep = 
i;F/(2Im£!j?(fc, £fc)), as given in Sec. fll] Thus, we 
first present the theoretical formalism for calculating 
Im£f p (fc,w). We consider the electron-phonon interac- 
tion through the deformation potential coupling. In this 
case, Im£g P (fc,cj) is given by-2^: 

s f d 2 Q , ,9 1 + ss' cos 9 
ImEf (*,«)= » £ / 7^I M I 2 

x [n F (u) q + uu)) + n B (u) q )] S(lj + vui q — e'), (23) 

where \M\ 2 = ® ^ is the electron-phonon scattering 
matrix element, while D, p m and vi are the deformation 
potential, the graphene mass density and the phonon ve- 
locity, respectively, s, s' — ±1 denote the band indices, 
and v = 1 ( — 1) corresponds to the absorption (emis- 
sion) of an acoustic phonon with frequency u> q = viq. 
e' = s'Kvf^ + q| — yU, and 6 is the angle between k and 
k + q. 



B. The asymptotic behavior of ImEs P (fc,^) 

We start with the asymptotic behavior for the 
energy- and temperature-dependence of Im£!p, cal- 
culated analytically within the quasielastic scattering 
approximatio n 30 ' 31 , namely |k| ~ |k + q|. This approx- 
imation is justified by the fact that the phonon velocity 
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is much smaller than the Fermi velocity vf in graphcnc. 

At zero temperature, np becomes the Heaviside unit 
step function, ns(cj q ) = 0, and the interband scattering 
vanishes. Then, the integration in Eq. (|23[) can be carried 
out in the limit <C Ep as: 



Im£m,&) 



1 



D 2 k F 



VF 

8n p m V[V F \ Vi 
ImE2 > (fc,^ fe 



E F 



(24) 



which shows that ImE_^(fc, ^) increases quadratically 
with increasing quasiparticle energy, in contrast to 
Im£ + (fc,£fc) oc £|ln£fc induced by e-e interaction. 
Eq. (HU) implies that l ep cx ^ 2 at T = 0. 

Next, we consider the temperature dependence of 
ImE^(fc, In the low temperature Bloch-Griineisen 
limit, T <C Tbgj with T^g = 2kpvi/ks, where the 
phonon system is degenerate, we get the following asymp- 
totic form for the on-shell ImS^(fc,w): 



10 cm t = 100 K 
10 10 cm^ 2 T =200 K: 
10" cm" 2 T= 100 K 
lo" cm" 2 T = 200K 




n(10 



tt n 2 k 2 



(—) <*) 

2 PmVlVF \T B G J 



In the Bloch-Griineisen regime, the temperature depen- 
dence of l ep is a power law with l ep cx T~ 2 . We also find 
that l ep does not depend on the carrier density in the low 
temperature limit. 

On the other hand, in the high temperature regime 
Tp ^> T ^> Tbg, where the phonon follows the nonde- 
generate equipartition distribution, the asymptotic for- 
mula for the on-shell ImE^f(fc,w) is given by: 



lmE^(kp,^ F ) 



1 D 2 k 2 F 



2 p m viVF T B g 



(26) 



Eq. (|26p shows that l ep decreases inverse linearly with in- 
creasing temperature, and is proportional to n~ x l 2 . We 
note that Tbg ~ ^f ~ \fn has a carrier density depen- 
dence and increases with increased doping in the system. 



C. Numerical results of l ep 

Fig. [8] shows our numerical results for the mean free 
path l ep . Fig. [Hlja) plots l ep as a function of tempera- 
ture for different values of carrier density. It is clear that 
l ep decreases monotonically with increasing temperature. 
Fig. [SJb) and (c) demonstrate that l ep is a decreasing 
function of and the carrier density n. From Fig. [SJc), 
we see that l ep is longer than 1000 nm for T < 100 K and 
n < 10 13 cm' 2 . Comparing Fig. [5] (a) with Fig. [5] (a), we 
find that l ep is about one order of magnitude longer than 
the inelastic mean free path / induced by e-e interaction 
for T < 100 K and n < 10 12 cm" 2 . Thus, the inelastic 
scattering is dominated by electron-electron scattering 
processes in this regime. With increasing carrier density 
and temperature, l ep becomes comparable to I, implying 
a considerable contribution of e-p interaction to the in- 
clastic scattering processes. We mention that in the pres- 



FIG. 8. (Color online). Calculated mean free path l ep induced 
by electron-phonon interaction. We use^i D — 25 eV, vi = 
2.6 x 10 6 cm/s and p m = 7.6 x 1(T 8 g/cm 2 . (a) l ep for £ k = 
as a function of temperature T for different carrier densities, 
(b) l ep as a function of energy ^/Ef for different carrier 
densities and temperatures, (c) l ep for £fc = £k. F as a function 
of carrier density for different temperatures. 



ence of both e-e and e-p interaction, the net quasiparticle 
scattering rate is given by l/r ee + 1/T ep and thus the net 
inelastic mean free path goes as ~ l/(l/l ee + l/l ep )- 



V. DISCUSSION AND CONCLUSION 

Before concluding, we first discuss the similarity and 
difference between graphene and conventional 2DES. Be- 
cause of the gapless nature of graphene, there are two 
contributing bands in graphene, i.e., conduction band 
and valence band, which allow both interband and intra- 
band single particle excitations. By contrast, the conven- 
tional 2DES has only one contributing band, either the 
conduction or the valence band, because of the large en- 
ergy gap. In addition, at zero temperature, ImE5(k, £k) 
of graphene is a smooth function of £kj where both plas- 
mon and interband excitations are absent. While there 
is a discontinuity in ImS R (k, £k) for a 2DES, caused by 
the plasmon^i and electron-hole excitation a 21 ' 32 , there is 
no such discontinuity in graphene due to its gaplessness. 

We now discuss the difference between extrinsic and 
intrinsic graphene. The low-energy quasiparticle scatter- 
ing rate of extrinsic graphene is quite similar to conven- 
tional 2DES, which is well described by the Fermi liquid 
theory. The non-monotonic temperature dependence of 
ImS^(k, u>) in extrinsic graphcnc arises from the compe- 
tition between plasmon broadening effects and thermally 
smeared single particle excitation boundaries. On the 
other hand, the quasiparticle scattering rate of intrinsic 
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graphene is linear in energy uj, a signature of the marginal 
Fermi liquid behavior. The characteristic feature of the 
non-monotonic ImS^(k, uj) for intrinsic graphene as a 
function of temperature is due to an interplay between 
plasmon enhancement and enlarged single particle exci- 
tation continua. We have also compared the inelastic 
mean free path with the elastic mean free path in or- 
der to assess the density and temperature range where 
the inelastic scattering length is decisively longer so that 
quantum interference is in principle allowed in graphene. 
When I < l e , quantum interference induced localization 
effects cannot manifest itself. In addition, motivated by a 
recent experiment on double-layer graphene, we have ob- 
tained analytic results for inelastic scattering in double- 
layer graphene where one layer acts to screen the other 
layer. 

Our numerical calculation shows that inelastic scatter- 
ing due to electron-phonon interaction is negligible for 
T < 100 K, and e-e interaction dominates in this tem- 
perature range. The electron-phonon scattering starts to 
play a role in extrinsic graphene when the temperature 
is above 200 K. In particular, the temperature depen- 
dence of l ep changes from T~ 2 to T _1 as the temperature 
increases from the low-T Bloch-Gruneisen regime to the 
high-T equipartition regime. On the other hand, the low- 
temperature asymptotic temperature dependence of the 
electron-electron interaction leads to I ~ (T 2 lnT) _1 for 
extrinsic graphene and I ~ T _1 for intrinsic graphene. 

We should note here that our calculation of the inelas- 
tic scattering length based on the imaginary part of the 
graphene self-energy is carried out in the so-called "bal- 
listic limit" where disorder effects arc neglected in the 
calculation of the inelastic mean free path itself. As such 
our calculated inelastic mean free path is an upper limit 
on the phase breaking length 1$ applicable to graphene 
quantum interference phenomena. It is well-known that 
disorder has qualitative and quantitative effects on the in- 
clastic mean free path^ 3 - and the phase breaking length, in 
general suppressing the mean free path substantially from 
its ballistic limit. In the diffusive limit, we anticipate the 
inelastic mean free path due to e-e interaction to have the 
asymptotic low-temperature behavior of / ~ (TlnT) -1 



instead of the I ~ (T 2 In T)" 1 ballistic behavior we find in 
Sec. ITU The phase breaking length 1$ is anticipated to be 
Itj, ~ T~ x in the low-temperature diffusive regime. Thus, 
the actual inelastic length due to electron-electron inter- 
action would be smaller than the ballistic limit results 
we obtain in the current work. We note that the dis- 
order effect on the electron-phonon interaction induced 
inelastic mean free path l ep derived in Sec. IIVI is likely 
to be small 3 - 3 -, and our ballistic limit results should apply 
equally well to the diffusive regime also. 

To conclude, we have obtained analytical asymptotic 
behavior of ImS^(k, £k) in graphene utilizing the GW 
approximation. We have also analyzed ImE^k, £k) for 
a double-layer graphene system. We have shown nu- 
merical results for the inelastic mean free path I and 
ImE fl (k, uS). We emphasize that finite temperature has 
strong effects on the inelastic scattering mean free path. 
We have discussed the difference between graphene and 
conventional 2DES. We have also provided results for 
the inelastic mean free path arising from electron-phonon 
interaction which becomes important only at very high 
temperatures since the typical electron-phonon dimen- 
sionless coupling constant (< 0.04 for carrier density n < 
10 13 cm -2 ) is much less than the typical dimension- 
less coupling constant r s (^ 0.4) in graphene. Our use 
of the GW approximation, which is the leading order 
Feynman-Dyson perturbative expansion in the dynami- 
cally screened Coulomb interaction, should be an excel- 
lent approximation for graphene (with quantitative pre- 
dictive power) because of the small values of r s (< 1) in 
graphere on substrates. Going beyond the GW approxi- 
mation is a formidable task which may not be necessary 
at all in understanding graphene inelastic processes. 
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